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Abstract 

We introduce a general mathematical model of symbiosis between different entities 
by taking into account the influence of each species on the carrying capacities of the 
others. The modeled entities can pertain to biological and ecological societies or to 
social, economic and financial societies. Our model includes three basic types: sym- 
biosis with direct mutual interactions, symbiosis with asymmetric interactions, and 
symbiosis without direct interactions. In all cases, we provide a complete classification 
of all admissible dynamical regimes. The proposed model of symbiosis turned out to be 
very rich, as it exhibits four qualitatively different regimes: convergence to stationary 
states, unbounded exponential growth, finite-time singularity, and finite-time death or 
extinction of species. 



PACS: 02.30.Hq, 87.10.Ed, 87.23.Cc, 87.23.Ge, 87.23.Kg, 89.65. Gh 



Keywords: Models of symbiosis, Nonlinear differential equations. Dynamics of symbiotic 
systems 



1 Introduction 



The term symbiosis describes close and usually long-term interactions between different bio- 
logical species. Symbiotic relationships are well known in biological and ecological societies. 
Numerous examples can be found in the books [1-5] . Arguably, the symbiosis that is the most 
important to us as humans is the one between the typical human body and the multitudes of 
commensal organisms representing members of five of the six kingdoms of life. Specifically, 
the microbiome, with more than 2000 bacterial species and some 10^^ microorganisms (about 
ten times the number of cells of the human body), is a key partner to the human immune 
system, and to the human metabolism as it participates in the synthesis of essential vitamins 
and amino acids, as well as in the degradation of otherwise indigestible plant material, and 
of certain drugs and pollutants in the guts [6]. Let us also mention the associations between 
plant roots and fungi, that are a feature of many terrestrial ecosystems, and without which 
the kingdom of plants would not exist as we know it. 

Generally speaking, symbiosis is not restricted to biological systems. Many relations in 
social and economic societies can also be interpreted as examples of symbiosis. For instance, 
the economic and intellectual levels of human societies can be considered as symbiotic to 
each other. The interconnections between basic and applied research are symbiotic as well 
[7-9]. Economic and financial relations between different firms or between firms and banks 
can be treated as symbiotic. The notion of symbiosis has a very wide applicability to various 
relations, whether in biology, ecology, economy, or finance. Therefore, when we shall refer to 
societies, we keep in mind different kinds of societies, including biological, human, economic, 
and others, whose interactions can be considered from the point of view of symbiotic relations. 

In biology, one distinguishes three main types of symbiosis: mutualism, commensalism, 
and parasitism [1-5]. Extending the notion of symbiosis to social and economic systems, we 
meet the same types of symbiosis in many different forms and structures. As in biology, there 
can occur relations between different more or less independent parts or species, illustrated by 
the relations between plant roots and fungi [10], or in economics, the relations between firms 
and banks, or between economics and arts. Social systems are also often characterized by 
the existence of complicated relations, e.g., corresponding to the relations between a general 
system and a particular subsystem. Examples are the relations between the economics of a 
country and a particular branch of the economics, or between culture and language. Biolog- 
ical counterparts are widespread, as illustrated by the endosymbiotic origin of mitochondria 
and plastids of eukaryotic cells, by some marine annelid worms lacking mouth, gut and anus, 
and which rely on multiple extracellular bacterial endosymbionts for their excretory system 
[11], and more generally between the microbiome network and the mammal body [12]. 

Despite these ubiquitous examples, to our knowledge, there is no general mathematical 
model describing such symbiotic relations, whether in biological or social sciences. The aim 
of the present paper is twofold. First, we extend the notion of symbiosis to social systems. 
Under the name "social" , we mean various kinds of social relations, including economical and 
financial. Second, we suggest a general mathematical model of symbiosis characterizing the 
overall dynamics of this process for various systems. We concentrate in the present paper on 
the general properties of the considered model, only briefiy mentioning possible applications. 

The paper proceeds as follows. Section 2 presents the general structure of the model, 
and Section 3, a classification of different types of symbiotic interactions. Sections 4, 5, and 
6 focuse on the structure of the model of symbiosis in the presence of mutual, respectively, 
asymmetric, and in the absence of direct interactions. Section 7 concludes. 



2 General Model of Symbiosis 



Interactions between different co-existing species are usually modeled by equations of the 
Lotka-Volterra type [13,14]. Such equations are adequate for predator- prey relations between 
different species, whether one studies the simple co-existence of two species [13,14], or more 
complicated cases corresponding to high-dimensional dynamical systems, such as multiply 
connected food webs [15,16], or the inter-relations between different cells and pathogens 
inside biological organisms [17,18]. However, the symbiotic relations are known to be princi- 
pally different from the predator-prey relations [1-5] . Hence, they require another mathemat- 
ical representation. To our knowledge, there is no general mathematical model describing 
symbiotic relations. 

Our basic suggestion is to describe symbiotic relations by the mutual influence of the 
co-existing species on their respective carrying capacities. Indeed, in the presence of sym- 
biotic relations, the species act on the livelihood of each other by creating resources that 
others exploit, by using resources created by others without feedbacks, or by destroying the 
livelihood in the case of parasitism. In mathematical language, the livelihood is nothing but 
the carrying capacity. 

Let us consider several species, or system parts, quantified by their population number 
Ni, with the index i enumerating the species. For instance, in biology, iVj can be the 
number of individuals of the i-th species or, in finance, this can be the amount of money 
or, in economics, this could be the quantity of some goods. We start with the logistic-type 
equations 

^^=7.JV.-^. (1) 

dt " Ki ' ^ ' 

for the variables Ni = Ni{t) as functions of time i > 0. Here 7^ is the birth rate of the 
biological species i or the growth rate in economic systems. The coefficient Cj characterizes 
the intensity of mutual competition between the agents of the i-th species. We will thus 
consider only the case where 

7^ > , a > . (2) 

The principal difference from the standard logistic equation is that the carrying capacity is 
here considered to be a function 



Ki^Ai + BiSi{{Ni}) (3) 

of the quantities A^j, in order to account for the co-existing symbiotic species. The first term 
Ai is the carrying capacity of the given surrounding livelihood. The second term characterizes 
the carrying capacity produced by other species. We thus capture the typical feature of 
symbiotic relations via the mutual influence of each species on the carrying capacities of 
others. The symbiotic coefficient Bi defines the intensity of producing, or destroying, the 
carrying capacity in the process of symbiotic relations. We refer to Bi as the production 
coefficient, when B^ is positive, or as the destruction coefficient, when B^ is negative. Thence, 
we shall consider that 

A > , Bi^ (-00, 00) . (4) 

Since the nature of the mutual interactions is characterized by the sign of the symbiotic 
coefficient B^, the symbiosis function 



Si{{Ni}) > , 



(5) 



describing the type of symbiotic relations, can be taken non-negative. The mathematical 
structure of Eq. (1) with (3), describing symbiosis, is principally different from equations of 
the Lotka-Volterra type and cannot be reduced to the latter. 

As we mentioned in the introduction, in biology and ecology, three main categories of 
symbiotic relations can be distinguished (mutualism, parasitism, and commensalism) , de- 
pending on whether the influence of one species on another is positive, negative, or neutral 
[1-5]. This classification is straightforwardly linked to our model. Below, we explain this for 
the case of two species, when i = 1,2, which will be treated in what follows. 

Mutualism implies the relations in which both species extract some benefit from their 
relationship. In our model, this is equivalent to 

i?i > , B2 > {mutualism) . 

Parasitism means that one of the species can benefit, while the other is harmed in the 
process. Generally, one can say that parasitism is the relation in which at least one of the 
species is harmed. This means that one of the following inequalities is valid: either 

Si > , < , 

or 

Si < , S2 > , 

or 

Bi < , B2 < (parasitism) . 

Commensalism is the relation in which one of the species benefits, while another is not 
affected. Hence, one of the following inequalities is satisfied: either 

Si > , S2 = , 

or 

Si = , S2 > (commensalism) . 

In addition to this classification, it is possible to distinguish different kinds of mutual in- 
teractions embodied in the form of the symbiosis function Si({Ni}) obeying inequality (5), 
which will be specified later. 

Before specifying this function Si({Ni}), it is convenient to introduce dimensionless units. 
For this purpose, we may choose two characteristic scales, Neff and Zeff, to serve as the 
measuring units for the species variables A^i and N2, respectively. We thus define the dimen- 
sionless species characteristics 

X = , z = ^ (6) 
and the dimensionless carrying capacities 

= n AT ' = n y ' 

Let us measure time in units of I/71 and introduce the ratio 



(8) 



Then, for two symbiotic species, Eqs. (1) take the form 

dx dz 
dt Hi ' dt 1/2 

These equations, controlhng the time evolution of the functions x = x{t) and z = z{t), must 
be complemented by the initial conditions 

x(0) = xo , ^(0) = zo . (10) 

By definition, the quantities Ni, measuring the amount of the corresponding species, are 
positive. Hence, we are interested only in the non- negative solutions of Eqs. (9): 

x{t) > , z{t) > . (11) 

We then need to specify the forms of the symbiosis functions Si{{Ni}) obeying inequality 
(5) and, respectively, the forms of the carrying capacities (3). The forms of functions (5) 
can vary, depending on the kind of symbiotic relations between the considered species. 



(9) 



3 Classification of Forms of Symbiotic Relations 

There can exist three main types of symbiotic relations. 

A. Influencing Livelihood Through Mutual Interactions 

A common case is when the symbiotic species influence the livelihood of each other by 
means of mutual interactions. Known biological examples arc plant trees and mushrooms 
and many other plants and fungi [1-5]. In the social sciences, this could be the relationship 
between the economic level of a country and the intellectual level of society. Or this can 
be the relation between two firms, producing different kinds of goods in close collaboration 
with each other. Another example is the relation between basic and applied research [7-9]. 
The simplest form of the carrying capacities, representing livelihoods that are influenced by 
mutual interactions, can be written as 

Ki^Ai + BiNiN2 , ^2 = ^2 + B2N2N1 . (12) 

The first terms Ai and A2 define the given carrying capacities, which in dimensionless units 
are 

The second terms characterize the produced carrying capacities, corresponding to the mutual 
infiuence of symbiotic species on the livelihood of each other. The dimensionless symbiotic 
coefficients are 

^^71^^ g^H^, (14) 

Gi G2 
The dimensionless carrying capacities (7) become 

yi = oi + bxz , y2 = 02 + gxz . (15) 



The given carrying capacities are taken positive, 



ai > , 02 > , (16) 

which means that each of the two species can survive in the absence of the other one. The 
interaction coefficients b and g can be of any sign, depending on whether the relations are 
beneficial or parasitic. 

The scahng units N^ff and Z^fj can be chosen arbitrarily. It is convenient to choose 
them as 

^eff - , Zeff - . (17j 

Then one has 

Oi = 02 = 1 , (18) 

and the symbiotic coefficients are 

6=A_2_1, (19) 

Therefore, the carrying capacities (15) become 

1/1 = 1 + hxz , 1/2 = 1 + . (20) 



B. Influencing Livelihood Through Asymmetric Interactions 

Another possibility is when the interactions of symbiotic species are not symmetric, 
such that the carrying capacity of one of them is influenced by mutual interactions, while 
the carrying capacity of another is influenced solely by the other species, not involving 
their mutual interactions. This type of symbiosis is common to many biological as well as 
social systems, when one of them is a subsystem of a larger one. For example, the relation 
between the economic level of a country and the development level of a particular branch of 
economics, say between the country gross domestic product and the level of science; or the 
relation between culture and language; or the relation between the size of social groups and 
the size of brain [19,20]. In that asymmetric case, the carrying capacities are 

/Ti = + EiAiAa , 7^2 = ^2 + ^2^1 . (21) 

In dimensionless units, a\ and 02 have the same form as in Eq 
coefficients are 

Ci ' Zeff 

The carrying capacities (7) are given by the expressions 

yi = oi + bxz , y2^a2 + gx . (23) 

Opting for the scaling units (17) yields normalization (18), with the symbiotic coefficients 

^^7^81^2 
C1C2 



(13), while the symbiotic 

(22) 



71^1^2 

C1A2 



(24) 



Thus, the dimensionless carrying capacities (23) become 



7/1 = 1 + bxz , 1/2 = I + gx . (25) 
C. Influencing Livelihood Without Direct Interactions 

This is probably the most common of biological and ecological symbiosis. Numerous 
examples can be found in the books [1-5]. In this case, the livehhood of each species is 
influenced by the presence of another species without involving their direct interactions. 
The carrying capacities are correspondingly given by 

K, = Ai + B1N2 , K2 = A2 + B2N1 . (26) 

The terms Ai and A2 are the a priori given carrying capacities, in the absence of the other 
species. Their dimensionless forms are the same as in Eq. (13), but the symbiotic coefficients 
are 

ClNeff C2Zeff 

The dimensionless carrying capacities (7) take the forms 

yi = ai + bz , y2 = a2 + gx . (28) 

Choosing again the same scaling units (17) gives the same normalization (18), but with the 
symbiotic coefficients 

''=^;cr- «=-c^- '''' 

This results in the dimensionless carrying capacities 

yi = l + bz , y2 = l + gx . (30) 

Note that here, as well as in the previous cases, the symbiotic coefficients b and g can 
take different values, depending on the strength of the mutual influence between the sym- 
biotic species. And these coefficients can be of different signs, describing cither beneficial 
or parasitic relations. In the following sections, we give a detailed analysis of each type of 
symbiosis classified above. 

4 Symbiosis With Mutual Interactions 

We study the equations (9), with the carrying capacities (20), leading to the coupled system 
of two ordinary differential equations: 

dx x^ dz z^ 

— = X — ; , — = Z — . (31) 

dt 1 + bxz dt 1 + gxz 

This system of equations is symmetric with respect to the change x ^ z and b ^ g. Three 
kinds of qualitatively different solutions are found for these equations, which we describe in 
turn. 

A. Convergence to Stationeu-y States 



Formally, there are five fixed points for Eqs. (31). However, as expressed by Eq. (11), 
we are looking for non-negative solutions. For each of the non-negative solutions, we ac- 
complish the Lyapunov stability analysis and select the stable solutions. Below, wc present 
the expressions and the stability analysis only for the stable stationary solutions that are 
non-negative. 

Equations (31) can possess a single stable fixed point given by the expressions 



1-6 + ^-^(1 + 6-5)2-46 ^ i + 5_^_^(i + ?,_^)2_45 
2g ' 26 ■ 

This point is stable in the stability region shown in Fig. 1 characterized by one of the 
conditions, either 

6 < , — oo < gi < oo , (33) 



or 



or 



0<6<1, 5<5c= (V6-l)'< 1, (34) 



6 > 1 , 5 < . (35) 
On the boundaries of the stability region, we have 

X* = 1 , = ^ (6 = 0, 5 < 1) , 

X* = ^ , z* = l (6 < 1, 5 = 0) , 

X* - ' " 71 (0 < 6 < 1, 5 = 5e) . (36) 

The solution to Eqs. (31) tends to the stationary point (32), provided that the symbiosis 
parameters 6 and g are in the stability region defined above, and when the initial conditions 
are in the attraction basin of this fixed point. For positive 6 and g, the attraction basin 
can be found only numerically while, if at least one of these parameters is negative, the 
attraction basin is characterized by one of the conditions, either 

x^zo < 1^ (6 < 0, 5 > 0) , (37) 

or 

xo^o < ^ (6 > 0, y < 0) , (38) 

or 

xqZq < min , (6 < 0, ^ < 0) . (39) 

The approach to the stationary solution (32) can be either monotonic or non-monotonic, from 
above or from below, depending on the parameters 6 and g and on the initial conditions xq 
and zq. In Figs. 2, 3, and 4, the typical behavior of the solutions is shown for different 
parameters and initial conditions: when xq > x* and zo > z* (Fig. 2), when xq > x* but 
Zq < z* (Fig. 3), and when xq < x* with zq < z* (Fig. 4). The case when xq < x* but 
Zq > z* is similar to the case when xq > x* but Zq < z*, with changing x hy z and 6 by g. 



It is instructive to compare the behavior of the symbiotic solutions, described by the 
coupled Eqs. (31), with that of the solutions of the decoupled equations {b — g — 0) 

^ = x-x\ ^^z-z^ (b = g = 0). (40) 

The solutions for these non-symbiotic species tend to the stable fixed point x* = z* = 1. 
The comparison, explicitly showing the role of symbiosis, is demonstrated in Figs. 5 and 
6. These two figures illustrate the general property that beneficial (respectively, parasitic) 
symbiosis leads to the increase (respectively, decrease) of the stationary solutions. Another 
important message is that stationary states can exist even in the presence of parasites. 

B. Unbounded Exponential Growth 

The solutions to the symbiotic equations (31) grow to infinity with increasing time t ^ oo, 
when there are no stable fixed points, that is, when either 

< 6 < 1 , g> Qc, (41) 

or 

6 > 1 , g>0. (42) 

A similar exponential divergence at infinity occurs when the stable fixed points exist, but 
the initial conditions are taken outside of the attraction basin, so that either 

xozo > ^ {b<0, g>0), (43) 

or 

XqZo > ^ {b>0, g<0) , (44) 
1^1 

or 

xqZo > max 1^ , ^| (b < 0, g < 0) . (45) 

The typical behavior of such increasing sohitions is shown in Fig. 7. When the initial condi- 
tions are outside of the attraction basin, this means that at least one of them is sufficiently 
large, that is, the initial population is so large that even the existence of parasites cannot 
suppress the development of this species. 

C. Finite-Time Death and Singularity 

A specific behavior occurs under parasitic symbiotic relations, when stable fixed points 
can exist, but the initial conditions are taken outside of the attraction basin, so that 

^ < xozo < ^ (6 < ^ < 0) . (46) 

Then, at the critical time tc, defined as the solution of the equation 

xitMQ = ^ , (47) 

the first variable x{t) sharply drops to zero, while the second one rises to infinity: 



x{t) , 



x{t) — >■ — oo 



{t ^ tc) , 



z{t) +00 , z{t) +00 {t tc) . (48) 



Here the overdot means time derivative. 

The behavior is inverted for the case when 



^ < XqZq < ^ (5 < 6 < 0) , (49) 

|o| \g\ 

for which the critical time is given as the solution of the equation 

x{Qz{tc) = ^ . (50) 
At this time, the first solution rises to infinity, while the second one quickly dies: 

x{t) +00 , x{t) +00 {t tc) , 

z{t) , z{t) -00 {t tc) . (51) 

The situation corresponding to the case (46) is illustrated in Fig. 8. The meaning of this 
phenomenon, when one of the species dies, while the other rises, is easily understandable. 
Under the mutually parasitic relations, each species destroys the livelihood of the other one. 
That species, whose hvelihood is destroyed faster, dies out, while the other species grows at 
the expense of the former one. 

5 Symbiosis With Asymmetric Interactions 

If the symbiotic relations between species arc characterized by asymmetric interactions with 
the carrying capacities (25), then Eqs. (9) take the form 

dx x'^ dz z^ . ^, 

— = X — , . (52) 

dt 1 + bxz dt 1 + gx ^ ^ 

The following qualitatively different dynamic regimes can occur. 

A. Convergence to Stationciry States 

There can exist just one stable fixed point, given by 



^ l-b-^{l-b)^-4bg ^ l + b-^{l-b)^-4bg 

X — —TT, ) z — . (53) 

2bg 2b ^ ^ 

Let us introduce the critical parameter value 

The fixed point (53) is stable and non-negative if either 

b<-l, g>gc, (55) 

or 

-1< & < , g>-l, (56) 



or 

< 6 < 1 , -l<9<gc, (57) 

or 

6 > 1 , -1 < c/ < . (58) 

The overall region of stability is depicted in Fig. 9. 

The marginal values of the fixed point, occurring at the boundary of the stability region, 

are 

= z*^l + g (6 = 0,^>-l), 

X* = ^ , ^* = 1 (6 < 1, ^ = 0) , (59) 

if 5r 7^ 5'c, and 

i9-9c), (60) 

ii g = and either 6 < — 1 or < 6 < 1. 

The solutions tend to the fixed point, provided that the related initial conditions are 
in the basin of attraction. When one of the symbiotic coefficients is negative, the basin of 
attraction is defined so that either 

xoz^ < 1^ (& < 0, 5 > 0) , (61) 

or 

xo < ^ (6 > 0, ^ < 0) , (62) 

m 

or ^ ^ 

Xo< — , XqZq <— (6 < 0, 5 < 0) . (63) 

The convergence to the stationary state can be monotonic or non-monotonic, as illustrated 
in Figs. 10 and 11 for various initial conditions and different symbiotic coefficients. The 
convergence of solutions to the stable fixed point is qualitatively similar to that documented 
in the case of symmetric interactions. 

B. Exponential Growth at Infinity 

Another regime that is similar to that found for symmetric interactions is the exponential 
growth at increasing time i — > oo. This happens when there are no stable fixed points and 
either 

< 6 < 1 , g>g,, (64) 

or 

6 > 1 , (? > , (65) 

with the critical parameter gc given in Eq. (54). The same kind of exponential growth arises 
when there exist stable fixed points, but the initial conditions are outside of their basin of 
attraction, which occurs for 

0<6<1, Q<g<gc, (66) 

or when 

xoZq > ^ (6 < 0, ^ > 0) . (67) 





The temporal behavior is quahtatively similar to that of Fig. 7. 



Finite-Time Singularity 

A different regime occurs when the parameter g is negative, and no stable fixed point 
exists, which occurs when 

0<6<oo, g <-l , (68) 

or 

-1< 6 < , g<-l, (69) 

or if 

& < -1 , g<gc<-l. (70) 
When the initial conditions obey the inequalities 

xo > , xozo > - J {b> 0, g <0) , (71) 

or the inequalities 

xo > ^ , xoZo > {b<0, g <0) , (72) 

then the solutions hve only a finite life until the critical time tc- The value of this critical 
point can only be determined numerically. As time approaches tc, the first variable remains 
finite, while the second one diverges: 

x{t) — >■ const > , x{t) ^ const {t — >■ tc) , 

z{t) oo , z{t) ^ oo (t^tc) . (73) 

The main difference with the case of Eq. (48) is the finiteness of the solution x{tc) at the 
critical time. The typical behavior of the two variables in this regime is shown in Fig. 12. 
The explosive divergence of one of the species announces a change of regime: beyond tc, the 
system of two species will change their characteristics and the nature of their interactions. 

The same kind of behavior also happens when the stable fixed points exist for b and 
negative g in the region 

b> -1 , -1 < g <0 , 

or 

6 < -1 , gc<g<0, 

but initial conditions are taken outside of the basin of attraction, so that Xq, Zq are defined 
by Eqs. (71) or (72), correspondingly. 

D. Finite-Time Death 

An interesting behavior appears for some parameters, in particular, when the parameter 
g is negative, so that no stable fixed points exist or when fixed points exist but the initial 
conditions are outside their domain of attraction. For g <Q, when the symbiosis coefficients 
are taken outside of the stability region presented in Fig. 9, and initial conditions obey the 
following inequalities, either 

xo < ^ , xqZq > - \ {b>0, g <0) , (74) 



or 

Xo < ^ , XqZq < (6 < 0, ^ < 0) , (75) 

or 

1 1 

Xo<-- , XqZq > — (b<0, g<0) , (76) 

or 

Xo> — , xqZq <— {b<0, g <0) , (77) 
1^1 \b\ 

then the solutions have a finite hfe and one of the species dies or exhibits a gradient catas- 
trophe at the death time ta (see below) found as the solution of the equation 

\gHtd) = 1 , (78) 

at which either z{t) = — oo or x{t) = — oo for t — >■ td- 

When the parameters are taken outside of the stability region and the initial conditions 
obey inequalities (74) or (75), then 

x(t) -— , x(t) const (t td) , 

m 

z{t) ccmst , z{t) -oo {t td) . (79) 
When the initial conditions obey inequality (76), then 

x{t) 1^ , x{t) +00 {t td) , 

z(t)->|^, i(t)^-oo {t^td), (80) 
while, when the initial conditions satisfy inequality (77), then 

x{t) 1^ , x{t) -oo {t td) , 

z{t) Ifl , z{t) ^ +00 {t^ td) . (81) 
\b\ 

In the cases (74) and (75), the solutions (79) remain non-zero at t — >■ td, but the time 

derivative of z{t) tends to — oo. Such a behavior corresponds to the so-called gradient 
catastrophe. The value z{td) is close to zero, as shown in Fig. 13. This kind of behavior 
can be interpreted as the extinction of the species described by the variable z{t). Therefore, 
after td, we can set 

z{t) = {t>td). (82) 

Then the continuation of the evolution of the population x{t) after td is characterized by 
the single equation x — x — x'^. The death of the species z{t) is due to the large negative 
symbiosis coefficient g, which leads the species x to suppress the population z of the other 
species, until its complete demise occurring at the time td- 

A behavior similar to that, described by Eqs. (80) and (81), occurs for g < 0, when the 
symbiosis coefficients are in the stability region, but the initial conditions obey inequalities 
(76) or (77). 



When {xo,zo} obey Eq. (76), then again there is a death time ta, given by the same 
Eq. (78). Here z{td) is not necessarily close to zero, but it experiences the same gradient 
catastrophe as in Eq. (79). Contrary to the case of Eq. (79), the derivative x{t) tends to 
+00. Strictly speaking, the coupled system of Eqs. (52) has no solution after the time ta- 
However, we can again interpret the gradient catastrophe, with z{t) — )■ — oo as the death of 
the species z{t), extending the solution x after ta using condition (82), as is shown in Fig. 
14. 

In the case when the initial conditions obey Eq. (77), the evolution also ends at the 
death time td, defined by the same Eq. (78). At this point, conditions (81) hold true. The 
difference with the case of Eqs. (80) is that here the gradient catastrophe happens with 
x{t) —oo while z(t) — )■ +oo for t — )■ td- Hence, interpreting td as the time of extinction 
of the species z that arc killed by the species x, the subsequent evolution of the dynamics 
after td is obtained by setting 

x{t) = {t>td). (83) 

The overall behavior is presented in Fig. 15. Comparing Figs. 14 and 15, we see that the 
dynamics essentially depends on the initial conditions. 



6 Symbiosis Without Direct Interactions 

When symbiosis is characterized by the influence of species on the livehhood of each other 
without direct mutual interactions, the corresponding carrying capacities are given by Eqs. 
(30). Then, the symbiotic Eqs. (9) take the form 

ditij dz z , \ 

-r = ^- --, ;— , -r = z - • (84) 

dt 1 + bz dt l + gx ^ ^ 

These equations are symmetric with respect to the replacement b ^ g and x ^ z. The 
different possible regimes exhibited by Eqs. (84) are described as follows. 

A. Convergence to Stationary States 

There can exist the single stable fixed point 

X* - z* - (85) 

^ " 1 - 65 ' ' ' 1-bg' ^^^^ 

which is positive and stable if either 

-1 < 6 < , g>-l (86) 



or 

b>0, 0<g<g,= ^. (87) 

The corresponding stability region is shown in Fig. 16. On the boundary of the stabihty 
region, the fixed point degenerates to one of the values 

X* = , z* = 1 {b = -1, g > -1) , 



X* = 1 , ^* = (6 > -1, ^ = -1) , 



X* , 0<z* <1 (b^g^-l) . (88) 

For 6 > and g > 0, the basin of attraction is the whole region in the b — g plane, where 
inequalities (87) arc valid. When one of the parameters 6 or is negative, or both are 
negative, then the attraction basin is defined by one of the following conditions, either 

a;o < ^ (6 > 0, 5 < 0) , (89) 

or 

^0 < ^ {b<0, g>0), (90) 

or ^ ^ 

Xo<-- , Zo<jTT {b<0 g<0) . (91) 

1^1 \b\ 

The convergence to the stationary solution can be monotonic or non-monotonic, depending 
on the symbiosis parameters and initial conditions, similarly to Figs. 2, 3, 4, and 10, 11. 

B. Exponential Growth at Infinity 

This behavior, which is analogous to that displayed in Fig. 7, occurs when 

b>Q, g>g,= ^. (92) 

Then, both solutions x{t) and z{t) exponentially increase as t ^ oo. 

C. Finite-Time Singularity 

When one or both of the symbiosis parameters are negative, implying parasitic symbiotic 
relations, a finite-time singularity can occur. For the initial conditions 

xo > ^ , -^0 > - ^ (6 > 0, ^ < 0) , (93) 

there exists a critical time tc, such that 

x{t) — >■ const , x{t) — >■ const {t — >■ tc) , 

z{t) +00 , z{t) +00 {t tc) . (94) 

This behavior is similar to that shown in Fig. 12. Because of the equation symmetry, there 
is the opposite case, occurring for 

xq> - -^0 > 1^ {b<0, g>0) , (95) 

when 

x{t) +00 , x{t) +00 {t tc) , 
z{t) — >■ const , z{t) — >■ const {t — >■ tc) ■ (96) 
If both symbiosis parameters are negative, and 

xo > ^ , -^0 > 1^ (6 < 0, 5 < 0) , (97) 



there exists a finite-time singularity of tlie same type as above, with the particular behavior 
depending on the parameter values and initial conditions. 



D. Finite-Time Death 

The occurrence of a gradient catastrophe, of the type described previously in Sec. 5, 
depends on the values of the symbiotic parameters and initial conditions. If the symbiotic 
parameters are outside of the stability region and the initial conditions are such that 

< ^ , ^o>-l (6 > 0, 5 < -1) , (98) 

then there exists a finite time td, determined as the solution of the equation 

\9\x{td) = 1 , (99) 

when ^ 

x{t) — )■ j— I , x{t) — >■ const {t — >■ td) , 

z{t) const , z{t) -oo {t td) ■ (100) 

The gradient catastrophe occurs for the species z, whose gradient tends to — oo. This can 
be interpreted as an abrupt collapse of 2; to and thus as the extinction of this species. The 
continuation of the solution x is then obtained by setting z = ior t > td- The extinction 
of the species z is caused by the strong parasitic action of the species x onto species z, 
represented by the negative value of the symbiotic parameter g. 
Under the condition 

Xo>--, zo<^ {b<-l, g>0), (101) 
9 \b\ 

the role of x and z is interchanged and the finite-time extinction is observed for the species 
X. At the finite time td, given as the solution of the equation 

\b\z(td) = 1 , (102) 

we have 

x{t) — >■ const , x{t) — >■ —00 {t td) 

z{t) ^ , z{t) canst {t td) . (103) 

Again, interpreting the gradient catastrophe for x as its extinction, the continuation of the 
solution for z is obtained by setting x = ior t > td- The extinction of the species x is 
caused by the strong parasitic action of z on x, due to the negative symbiotic parameter b. 

The gradient catastrophe can occur even when the symbiotic parameters are inside the 
stability region, but with initial conditions that are outside of the basin of attraction. Thus, 
if 

1 1 

Xo<^ , Zo> Tji {b<0, g<0) , (104) 
then there exists a finite time td, given as the solution of the equations 

\g\x(td) = \b\z(td) = 1 , (105) 



at which ^ 

x{t) — >■ j— I , i; — >■ +00 {t — >■ td) , 

z{t) ^ ^ , (t^td) . (106) 

Again, here the gradient catastrophe for z can be understood as its extinction. Conversely, 
when the initial conditions are such that 



xo>r^, zo<jjj {b<0, g<0), (107) 



1 


1 


1 r ' 


^0 < TIT 








then there is the finite time td, given by the same Eqs. (105), when 

m 

^ii) ^ ' ^(t) ^ +00 {t^ td) . (108) 
\b\ 

This can be interpreted as the extinction of the species x, after which it can be set zero for 
t > td. 

The overall behavior is presented in Fig. 17. 

7 Conclusion 

We have proposed a general approach to describe the dynamics of symbiotic relations between 
several species. We have argued that the notion of symbiosis can be apphed not merely to 
biological systems but can also be generalized to social systems of different nature. 

The principal point of our approach is the description of symbiosis through the influence 
of each species on the carrying capacities of the others. This is in agreement with the common 
understanding that symbiotic species act on the livelihood of each other, either improving 
it, under beneficial relations, or destroying it, under parasitic relations. 

The general symbiotic model can be represented by several variants characterizing dif- 
ferent types of symbiotic relations. We have considered three basic types of such relations, 
symbiosis with mutual interactions, symbiosis with asymmetric interactions, and symbiosis 
without direct interactions. In all cases, we have provided a complete classification of all 
admissible dynamical regimes. 

The functional dependence of the carrying capacities on the species variables x and z 
constitutes the principal difference between our model and the logistic equation, where the 
carrying capacity is fixed. We have chosen here the simplest functional form of the carrying 
capacity characterizing symbiotic interactions. The carrying capacity has been taken as a 
combination of terms not exceeding the bilinear order in terms of the species variables. In 
general, it would be possible to describe the carrying capacity as an expansion in increasing 
powers of x and z. Then, we could use an effective summation of the expansion, obtaining 
a more complicated expression, e.g., in the form of the self-similar exponentials, as has 
been done for the model describing the dynamics of a nonequilibrium financial system [21]. 
Another modification could be by including the effects of delays into the carrying capacity, 
as has been done for the model of punctuated evolution [22] . 



Even in its simplest form, the proposed model of symbiosis turned out to be sufficiently 
rich, exhibiting four qualitatively different regimes: convergence to stationary states, un- 
bounded exponential growth, finite-time singularity, and finite-time death or extinction of 
species. The characteristics of real systems, of course, cannot diverge. Therefore, the occur- 
rence of finite-time singularities should be understood as the manifestation of an instability 
developing in the system, which signals a change of regime into a new structural phase 
described by different mechanisms and thus different equations [23]. In some cases, the sin- 
gularity can be avoided by slightly modifying the model, c.g, by taking into account higher 
order powers of the variables x and z. Such higher orders can remove the finite-time singu- 
larity [21]. In any case, when a finite-time singularity does happen, this can be understood 
as the signal that the system behavior is drastically changing at this point, somewhat similar 
to the changes occurring under phase transitions [24,25]. 

The suggested model of symbiosis can be applied to a variety of systems, biological, social, 
economic, financial, and so on. We have just mentioned some of the possible applications. 
The main goal of the present paper has been to advance a general parsimonious model 
of symbiosis and to analyze its main dynamical regimes. Here, we limit ourselves by the 
mathematical side of the problem. Particular applications require separate investigations 
and will be studied in future publications. 
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Figure Captions 



Fig. 1. Region of stability (shaded) in the parameter plane b — g for the fixed points in 
the case of symbiosis with mutual interactions. 

Fig. 2. Convergence to stationary states of the solutions x{t) (solid line) and z{t) 
(dashed-dotted line), as functions of time, in the case of symmetric mutual interactions, 
when xq > x* and zq > z*, for different parameters b and g and initial conditions {xo,zo}: 
(a) b = -0.1, 5 = 1, {0.75,3.8}i, {1,3.8}2, {2.5,3.8)3; the fixed points being x* = 0.730, 
z* = 3.702. (b) b = 0.5, g = 0.01 < gc = 0.086, {2.5, 1.5}i, {4, 1.5)2, {7, 1.5)3; x* = 2.043, 
z* = 1.021. (c) b = -0.5, g = -0.1, {0.8, l)i, {0.8,1.5)2, {0.8,2.499)3, x* = 0.681, 
z* = 0.936. (d) b = 0.5, g = -0.1 <gc = 0.086, {3, l)i, {3, 3)2, x* = 1.742, z* = 0.852. 



Fig. 3. Convergence to stationary states of the solutions x{t) (solid line) and z{t) 
(dashed-dotted line), as functions of time, in the case of symmetric mutual interactions, 
when Xq > X* but Zq < z*, for different parameters b and g and initial conditions {^c^^o)- 
(a) b = 0.05, g = 0.5 < gc ^ 0.603, {4,0.05)i, {4,0.5)2, {4,2)3; the fixed points being 
X* = 1.130, z* = 2.298. (b) b = 0.7, g = 0.02, {6,0.05}i, {6,0.5)2, {6,0.9)3; x* = 4.258, 
z* = 1.093. (c) b = -l,g = 0.1, {0.45, 0.01)i, {1, 0.01)2, {1.5, 0.01)3, x* = 0.438, z* = 1.281. 
(d) b = -0.01, g = l, {3, 0.05)1, {3, 2)2, {3, 3.33)3, X* = 0.730, z* = 3.702. 



Fig. 4. Convergence to stationary states of the solutions x{t) (solid hne) and z{t) 
(dashed-dotted line), as functions of time, in the case of symmetric mutual interactions, 
when Xq < X* and zq < z*, for different parameters b and g and initial conditions {xq,zq}: 
(a) b = 0.1, g = 0.45 < gc = 0.468, {0.05, 2)i, {0.5,2)2, {1.1,2)3; the fixed points being 
X* = 1.333, z* = 2.5. (b) b = 0.7, g = 0.025 < = 0.027, {0.1,0.01)i, {1.4,0.01)2, 
{4.9,0.01)3; X* = 5, = 1.143. (c) 6 = -1, 5 = 0.1, {0.01,0.1)1, {0.2,0.1)2, {0.487,0.1)3, 
X* = 0.4874, z* = 1.051. (d) b = -0.1, g = -1, {0.01, 0.001)i, {0.01,0.3)2, {0.01,0.51)3, 
X* = 0.951, z* - 0.513. 

Fig. 5. Comparison of the symbiotic solutions x{t) (solid line) and z{t) (dashed-dotted 
line) with the solutions x{t) — z{t) (dashed line) of the decoupled equations (40) for the 
same initial conditions = ^0 = 0.1 < 1, but for different symbiotic parameters b and g: 
(a) b = 0.25, g = 0.1 < gc = 0.25; the stationary points of the symbiotic equations being 
X* — 1.411, z* — 1.164. (h) b = 2, g — —0.5, the fixed points of the symbiotic equations 
being x* — 3.562, z* — 0.360. (c) b — —1, g — 2, the symbiotic fixed points being x* — 0.293, 
z* = 2.414. (d) 6 = -1, 5 = -2, the symbiotic fixed points being x* = 0.707, z* = 0.414. 

Fig. 6. Comparison of the symbiotic solutions x{t) (solid line) and z{t) (dashed-dotted 
hne) with the solutions x{t) — z{t) (dashed hne) of the decoupled equations (40) for the 
same initial conditions = ^0 = 1-8 > 1, but for different symbiotic parameters b and 
g: (a) b — 0.5, g — 0.080 < gc — 0.086; the fixed points of the symbiotic equations being 



X* = 2.823, z* = 1.292. (b) b = 0.5, g 
being x* = 1.742, z* = 0.852. (c) b 
X* = 0.582, z* = 2.393. (d) b = -0.1, g 
z* = 0.782. 



= —0.1, the fixed points of the symbiotic equations 
= —0.3, g = 1, the symbiotic fixed points being 
= —0.3, the symbiotic fixed points being x* = 0.927, 



Fig. 7. Logarithmic behavior of the exponentially growing solutions x{t) (solid line) 
and z{t) (dashed-dotted line) for different symbiotic parameters and initial conditions: (a) 
h = l^g = 0.1, xq = 148, Zq = 0.135. (b) b = 0.5, g = l> g^ = 0.086, Xq = 0.368, Zq = 0.05. 
(c) 6 = 1, ^ = 0.1 > g^^O,xo^ 0.135, zo = 2.72. (d) b ^ -1, g ^ 0.1, Xq = 148, Zq = 0.05. 

Fig. 8. Logarithmic behavior of the solutions in the case of finite-time death and 
singularity, x{t) (solid line) and z{t) (dashed-dotted line), for the parasitic relations with the 
symbiotic coefficients b — —1, g — —2, under the initial conditions Xq — 1.8, Zq — 0.5. For 
these parameters, the critical time is tc — 0.87. 

Fig. 9. Region of stability (shaded) in the parameter plane b — g ioi the fixed points in 
the case of symbiosis with asymmetric interactions. 

Fig. 10. Convergence to stationary states of the solutions x{t) (solid line) and z{t) 
(dashed-dotted line), as functions of time, in the case of asymmetric interactions, for the same 
initial conditions, such that Xq > Zq, with xq — 1.8 and Zq — 0.1, for different parameters 

b and g: (a) b = 0.5, g = 0.1 < gc = 0.125; the fixed points being x* = 2.764, z* = 1.276. 
(h) b = 2, g = -0.5; the fixed points being x* = 1.618, z* = 0.191. (c) b = -0.25, g = 2; 
with the fixed points x* = 0.638, z* = 2.275. (d) b = —0.3, g = —0.4; the fixed points being 
X* = 0.833, z* = 0.667. 



Fig. 11. Convergence to stationary states of the solutions x{t) (solid line) and z(t) 
(dashed-dotted line), as functions of time, in the case of asymmetric interactions, for the same 
initial conditions, such that xq < zq, with xq = 0.1 and zq — 1.8, for different parameters 
b and g: (a) b = 0.5, ^ = 0.1 < = 0.125; the fixed points being x* = 2.764, z* = 1.276. 
(h) b = 2, g = -0.5; the fixed points being x* = 1.618, z* = 0.191. (c) b = -0.25, g = 2; 
with the fixed points x* — 0.638, z* — 2.275. (d) b — —0.3, g — —0.4; the fixed points being 
X* = 0.833, z* = 0.667. 

Fig. 12. Logarithmic behavior of the solutions in the case of asymmetric interactions, 
in the presence of the finite-time singularity, x(t) (solid line) and z(t) (dashed-dotted line), 
for the symbiotic coefficients b = 2, g = —1.5 and the initial conditions Xq = 1.8, Zq = 0.1. 
The critical time here is tc = 4.43. 



Fig. 13. Finite-time death in the case of asymmetric interactions. The species x{t) (solid 



line) kill the species z{t) (dashed-dotted line) at the death time td- (a) b = —2, g = —1.5; 
the initial conditions are xq = zq = 0.1; the death time is td = 3.8. {h) b — 1, g — —1.5; the 
initial conditions are Xq = 0.1, Zq = 1.8; the death time is td = 2.6. 

Fig. 14. Finite-time death in the case of asymmetric interactions. The species x{t) 
(solid line) kill the species z(t) (dashed-dotted line) at the death time td- The symbiotic 
coefficients are b = —2, g = —1.5; the initial conditions are Xq = 0.1, Zq = 6; the death time 
is td = 0.559. 

Fig. 15. Finite-time death in the case of asymmetric interactions. The species z{t) 
(dashed-dotted line) kill the species x{t) (solid line) at the death time td- The symbiotic 
coefficients are b — —2, g — —1.5; the initial conditions are xq — 1.8, zq = 0.1; the death 
time is td — 1.016. 

Fig. 16. Region of stability (shaded) in the parameter plane b — g for the fixed points 
in the case of symbiosis without direct interactions. 

Fig. 17. Finite-time death in the case of symbiosis without direct interactions. Temporal 
behavior of solutions x{t) (solid line) and z{t) (dashed-dotted line) for different symbiosis 
parameters and initial conditions: (a) b — —0.75, g = —0.5, Xq = 0.8, Zq = 3; the death time 
being td = 0.204. (b) b = -1.5, g ^ 1, xq ^ 1, Zq ^ 0.1; with the death time td = 2.412. (c) 
b^ -1, g = -2, xo = 0.8, = 0.1; with the death time td = 4.279. (d) 6 = -2, ^ = -1, 

= 0.1, ^0 = 1; the death time being td — 1.459. 
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Figure 1: Region of stability (shaded) in the parameter plane b — g for the fixed points in 
the case of symbiosis with mutual interactions. 
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Figure 2: Convergence to stationary states of the solutions x{t) (solid line) and z{t) (dashed- 
dotted line), as functions of time, in the case of symmetric mutual interactions, when xq > x* 
and zq > z*, for different parameters b and g and initial conditions {xo,zo}: (a) b = —0.1, 
g = l, {0.75,3.8}i, {1,3.8}2, {2.5,3.8)3; the fixed points being x* = 0.730, z* = 3.702. (b) 
b = 0.5, g = 0.01 < gc = 0.086, {2.5, 1.5}i, {4, 1.5)2, {7, 1.5)3; x* = 2.043, z* = 1.021. (c) 
b = -0.5, ^ = -0.1, {0.8, 1)1, {0.8, 1.5)2, {0.8, 2.499)3, X* = 0.681, z* = 0.936. (d) b = 0.5, 
g = -O.K g^ = 0.086, {3, l}i, {3, 3)2, x* = 1.742, z* = 0.852. 







(a) 

z(t)^ 




/ 




\ 






x(t) 









2 






(b) 




















(d) 






z{t) 






x(t) 


■ (3) 




■ — 1 



Figure 3: Convergence to stationary states of the solutions (solid line) and z{f) (dashed- 
dotted line), as functions of time, in the case of symmetric mutual interactions, when xq > x* 
but Zo < for different parameters h and g and initial conditions {xq,z^: (a) 6 = 0.05, 
^ = 0.5 < = 0.603, {4,0.05}i, {4,0.5}2, {4, 2)3; the fixed points being x* = 1.130, 
z* = 2.298. (b) b = 0.7, g = 0.02, {6,0.05}i, {6,0.5}2, {6,0.9)3; x* = 4.258, z* = 1.093. 
(c) 6 = -1, ^ = 0.1, {0.45,0.01}i, {1,0.01)2, {1.5,0.01)3, X* = 0.438, z* = 1.281. (d) 
b = -0.01, c/ = 1, {3, 0.05)1, {3, 2)2, {3, 3.33)3, X* = 0.730, z* = 3.702. 




Figure 4: Convergence to stationary states of the solutions x(t) (solid line) and z(t) (dashed- 
dotted line), as functions of time, in the case of symmetric mutual interactions, when xq < x* 
and Zo < z*, for different parameters b and g and initial conditions {xo,zq}: (a) b = 0.1, 
g = 0.45 < gc = 0.468, {0.05, 2}i, {0.5, 2)2, {1.1,2)3; the fixed points being x* = 1.333, 
z* = 2.5. (b) b = 0.7, g = 0.025 < gc = 0.027, {0.1,0.01}i, {1.4,0.01)2, {4.9,0.01)3; 
X* = 5, z* = 1.143. (c) b = -l, g = 0.1, {0.01,0.1)i, {0.2,0.1)2, {0.487,0.1)3, x* = 0.4874, 
z* = 1.051. (d) b = -0.1, g = -1, {0.01, 0.001)i, {0.01,0.3)2, {0.01,0.51)3, x* = 0.951, 
z* = 0.513. 
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Figure 5: Comparison of the symbiotic solutions x(t) (solid line) and z(t) (daslied-dotted 
line) with the solutions x(t) = z(t) (dashed line) of the decoupled equations (40) for the 
same initial conditions xq = 2^0 = 0.1 < 1, but for different symbiotic parameters b and g: 
(a) b = 0.25, g = 0.1 < gc = 0.25; the stationary points of the symbiotic equations being 
X* = 1.411, z* = 1.164. (h) b = 2, g = —0.5, the fixed points of the symbiotic equations 
being x* = 3.562, z* = 0.360. (c) b = —1, g = 2, the symbiotic fixed points being x* = 0.293, 
z* = 2.414. (d) b= -1, g = -2, the symbiotic fixed points being x* = 0.707, z* = 0.414. 




Figure 6: Comparison of the symbiotic solutions x{t) (solid line) and z{t) (dashed-dotted 
line) with the solutions x{t) = z{t) (dashed line) of the decoupled equations (40) for the 
same initial conditions xq = -^o = 1-8 > 1, but for different symbiotic parameters b and 
g: (a) b = 0.5, g = 0.080 < gc = 0.086; the fixed points of the symbiotic equations being 
X* = 2.823, z* = 1.292. (b) b = 0.5, g = —0.1, the fixed points of the symbiotic equations 
being x* = 1.742, z* = 0.852. (c) b = —0.3, g = 1, the symbiotic fixed points being 
X* = 0.582, z* = 2.393. (d) b = -0.1, g = -0.3, the symbiotic fixed points being x* = 0.927, 
z* = 0.782. 



Figure 7: Logarithmic behavior of the exponentially growing solutions x{t) (solid line) and 
z{t) (dashed-dotted line) for different symbiotic parameters and initial conditions: (a) 6 = 1, 
g = 0.1, xo = 148, zo = 0.135. (b) b = 0.5, g = 1 > = 0.086, xq = 0.368, zq = 0.05. (c) 
b = l, g = 0.l> gc = 0, xo = 0.135, zq = 2.72. (d) b = -1, g = 0.1, xq = 148, zq = 0.05. 




Figure 8: Logarithmic behavior of the solutions in the case of finite-time death and sin- 
gularity, x{t) (solid line) and z(t) (dashed-dotted line), for the parasitic relations with the 
symbiotic coefficients b = —1, g = —2, under the initial conditions xq = 1.8, zq = 0.5. For 
these parameters, the critical time is tc = 0.87. 




Figure 9: Region of stability (shaded) in the parameter plane b — 
the case of symbiosis with asymmetric interactions. 
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Figure 10: Convergence to stationary states of the solutions x{t) (solid line) and z{t) (dashed- 
dotted line), as functions of time, in the case of asymmetric interactions, for the same initial 
conditions, such that xq > zq, with xq = 1.8 and zq = 0.1, for different parameters b and g: 
(a) b = 0.5, g = 0.1 < g^ = 0.125; the fixed points being x* = 2.764, z* = 1.276. (b) 6 = 2, 
g = —0.5; the fixed points being x* = 1.618, z* = 0.191. (c) b = —0.25, g = 2; with the fixed 
points X* = 0.638, z* = 2.275. (d) b = -0.3, g = -0.4; the fixed points being x* = 0.833, 
z* = 0.667. 
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Figure 11: Convergence to stationary states of the solutions x{t) (solid line) and z{t) (dashed- 
dotted line), as functions of time, in the case of asymmetric interactions, for the same initial 
conditions, such that xq < zq, with xq = 0.1 and zq = 1.8, for different parameters b and g: 
(a) b = 0.5, g = 0.1 < g^ = 0.125; the fixed points being x* = 2.764, z* = 1.276. (b) 6 = 2, 
g = —0.5; the fixed points being x* = 1.618, z* = 0.191. (c) b = —0.25, g = 2; with the fixed 
points X* = 0.638, z* = 2.275. (d) b = -0.3, g = -0.4; the fixed points being x* = 0.833, 
z* = 0.667. 
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Figure 12: Logarithmic behavior of the solutions in the case of asymmetric interactions, in 
the presence of the finite-time singularity, x(t) (solid line) and z{t) (dashed-dotted line), for 
the symbiotic coefficients b = 2, g = —1.5 and the initial conditions xq = 1.8, zo = 0.1. The 
critical time here is tc = 4.43. 




Figure 13: Finite-time death in the case of asymmetric interactions. The species x{t) (sohd 
hne) kill the species z(t) (dashed-dotted line) at the death time t^: (a) b = —2, g = —1.5; 
the initial conditions are Xq = zq = 0.1; the death time is td = 3.8. (h) b = 1, g = —1.5; the 
initial conditions are Xq = 0.1, Zq = 1.8; the death time is ta = 2.6. 




Figure 14: Finite-time death in the case of asymmetric interactions. The species x{t) (solid 
line) kill the species z{t) (dashed-dotted line) at the death time t^. The symbiotic coefficients 
are b = —2, g = —1.5; the initial conditions are Xq = 0.1, zq = 6; the death time is td = 0.559. 




Figure 15: Finite-time death in the case of asymmetric interactions. The species z{t) (dashed- 
dotted hne) kill the species x{t) (solid line) at the death time t^. The symbiotic coefficients 
are b = —2, g = —1.5; the initial conditions are Xq = 1.8, Zq = 0.1; the death time is 
td = 1.016. 




Figure 16: Region of stability (shaded) in the parameter plane b — g for the fixed points in 
the case of symbiosis without direct interactions. 




Figure 17: Finite-time death in the case of symbiosis without direct interactions. Temporal 
behavior of solutions x{t) (solid line) and z{t) (dashed-dotted line) for different symbiosis 
parameters and initial conditions: (a) b = —0.75, g = —0.5, xq = 0.8, zq = 3; the death time 
being ta = 0.204. (b) b = -1.5, g = 1, xq = 1, zq = 0.1; with the death time td = 2.412. (c) 
b = -1, g = -2, xo = 0.8, Zq = 0.1; with the death time td = 4.279. (d) b = -2, g = -1, 
Xq = 0.1, Zq = 1; the death time being td = 1.459. 



